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We present a family of generalized unitary transformations that simplifies the Hamiltonian for 
a harmonically trapped two level atom (or ion) interacting with a plane wave laser field. Novel 

(N . near resonant single as well as double vibrational phonon dynamical regimes are found. The valid- 

' ity condition of the often used motional rotating wave approximation (MRWA) is examined both 

numerically and analytically. Large errors are found within typical regimes of MRWA with respect 
to the motional degrees of freedom. The effects of MRWA in trapped ion systems are shown to be 
opposite to that of the rotating wave approximation (RWA) in the usual Jaynes-Cummings model. 
Our study points to a more restrictive condition on particle localization (Lamb-Dicke) parameter for 
the validity of MRWA in the single phonon dynamical regime. It also sheds new light on quantum 

l information storage and processing with trapped atoms. 
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In the last few years, much attention has been focused on quantum dynamics and coherence properties of trapped 
atoms or ions P|-^|. These studies led to many potentially attractive applications of harmonically trapped parti- 
cles, e.g. in generating nonclassical vibrational phonons Q, implementing fast quantum gates and in achieving 
phonon-ion entanglement J5J]. In most of these studies, a coherent plane wave laser field near-resonantly couples two 
electronic states of an atom. The inclusion of motional degrees of freedom leads to tremendous complication. Several 
simplifications have been developed that reduce a trapped particle dynamics to the familiar Jaynes-Cummings Model 
(JCM) type form with the use of Lamb-Dicke limit (LDL), the strong confinement limit, or the motional rotating 
wave approximation (MRWA). In fact, the JCM description can be achieved regardless of laser field configuration, i.e. 
whether it is a travelling- wave || or a standing- wave [0J^]. LDL requires the particle localization size, given by the 
harmonic trap ground state width a, to be much less than the near resonant laser wavelength A, i.e. r\ = (27r)a/A <C 1. 
It is often believed that MRWA for harmonically trapped particles works well within LDL as it demands a less 
restrictive condition 77 <§; 2 || or 77 <§; 4 ||. 
£NJ , In our ongoing effort to understand motional effects of trapped particles for quantum information processing [Q, 
this issue of validity regimes for trapped particle MRWA again arises. We note some earlier investigations show that 
rotating wave approximation (RWA) does not work satisfactorily for Hamiltonians containing multiple transitions [^|^] 
or for systems initially prepared in a superposition of internal states [fTofl . Furthermore, multi-particle properties such 
as entanglement can have different sensitivity dependence on single particle MRWA |llj . With harmonically trapped 
particles, this issue becomes particularly acute as the equal distant motional states span an infinite dimensional Hilbert 
space. MRWA is typically made after a linearization by either taking LDL or by applying an exact simplifying unitary 
transformation. Within the latter approach, errors from MRWA can be further modified by the back transformation 
'-^J when dynamical observables are calculated in the original frame. 

In this article we study quantum dynamics of a single trapped atom by employing a numerical diagonalization 
procedure without employing MRWA. We assess the validity regime for MRWA by comparing with results from 
analytical models obtained under MRWA. This paper is organized as follows. In Sec. II, our model system of a single 
trapped two level atom is described and a review of the unitary transformation method for linearizing the Hamiltonian 
is presented. The new result, the existence of a family of general unitary transformations for linearizing our model 
Hamiltonian, is then introduced. In Sec. Ill, we discuss two linearized models obtainable from our transformations. 
■ In Section IV, we outline several technical points about analytical solutions of the transformed model Hamiltonian 
under MRWA. We also discuss a numerical diagonalization procedure used for exact dynamic solutions. Selected 
results and comparisons are presented in Sec. V. Finally we conclude in Sec. VI. 

II. THE MODEL SYSTEM AND A FAMILY OF GENERAL UNITARY TRANSFORMATIONS 

To simplify our discussion, we consider a one-dimensional model of a harmonically trapped two level particle 
interacting with a near resonant laser field ||,[l2| . The system Hamiltonian is given by 
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where 0, u)l, and fci denote respectively the Rabi frequency, carrier frequency, and wave number of a coherent driving 
laser. The electron transition frequency between excited (|e)) and ground state (\g)) is uj eg . <r a b = \ a )(b\ (a/b = e,g) 
are atomic projection operators. Consist with convention we denote cr_ = \g)(e\ and a z = a ee — a gg . For a neutral 
particle, typically the approximate harmonic trap potential is internal state dependent, i.e. Vt a — (\/2)Mv 2 x 2 with 
a corresponding trap frequency v a . P is the motional momenta and M is the mass of the particle. 

We first simplify Eq. (m) by changing into the interaction picture (rotating frame) with the unitary transformation 
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where S 



ljl is the laser field detuning. We now introduce motional phonon annihilation (creation) operator 



for the ground state a (a^) according to P — y/%Mv g /2p and x — yjTi/2Mv g q with p 
The displacement operator becomes D(f3) = e /3a1_ ' 3 a , and Eq. (|J) simplifies to (fi = 1) 

$ a 2 v q + 5 
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where n = a^a is the phonon number operator and £ = (v 2 — v^)/Av g . The Lamb-Dicke parameter is now r\ — 
kh^/\j2Mv g . For this particular form of the Hamiltonian, we note that its interaction part (Hi) can be diagonalized 
by a general transformation matrix T = E^T, with 
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We choose fi £ 5i without loss of generality. It is easy to check that THiT 1 = (fi/2)a z , is independent of the arbitrary 
unitary functional operators Ei^ of a and a) . For any operator O, the transformed operator will be denoted to be 
= TQT^ . I — a ee + <jgg is the identity operator. This transformation can be compared to the generalized Power- 
Zienau transformation discussed earlier Jl3[ ] . In the phonon Fock state basis, it generates coherent superpositions of 
motional wave-packet states which was previously used in studying motional decoherence of atomic qubit operations 
JjJ. It reduces to a simplifying transformation used by Moya-Cessa et al. || when we take E\ = E% = D(—ir)/2). The 
T transformation on the Ho term can be conveniently calculated using the following properties 



D{a)D{p) = e (Q/3 *- Q * /3)/2 D(a + 0), 
D(^irf)a) aD(±ir]) = a) a ± irj((V — a) + J] 2 , 
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for arbitrary complex numbers a and (3. We note that Eq. (|J) simplifies to D(a)D(f3) = D(a + (3) when a and f3 
are purely imaginary. The second and the last term in the rhs of Eq. (^) correspond to the Doppler and recoil shifts. 
After this general transformation, we obtain 
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with g — r\Vgj2 and a x = a + + cr_. Different simplifications can be pursued by exploiting other forms of Ei 2 - We 
note that the first term in Eq. (^) is in the form of a squeezed displaced harmonic oscillator. Therefore it can be 
diagonalized in the Fock state basis by the squeezed coherent state transformation E% = E\ = D((3)S(^). To eliminate 
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the p dependence in the coefficient of I we choose (3 
p — > p — rj, and q — > q. We then obtain 



-ir]/2, which transforms according to n — > n — (j]/2)p + ry 2 /4, 
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The squeezing transformation S(£ = re lS ) causes a — > aa' — v*a with u = coshr and u = e l61 sinhr. For our model, £ 
is required to be real, i.e. 9 — 0. The action of S results in 
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Substitution these results into Eq. (||) leads to the elimination of the q 2 term in the coefficient of / term provided we 
choose r = (lne)/4 with e = 1 + 2Q/v g . This yields 
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We can now redefine parameters according to v 
transformed Hamiltonian 
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where the constant term in the coefficient of / has been dropped. We note that for ions with internal state independent 
trap frequencies £ = r = £ = 0. We recover the same result as in Ref. || by choosing (3 — —irj/2. The effects of different 
trap frequencies for different internal levels are a renormalization of energy parameters in the single vibrational phonon 
Hamiltonian and a double vibrational phonon interaction through a term quadratic in the position operator. Double 
phonon transitions are the typical interactions that can lead to squeezed and entangled states for the vibrational 
phonons. In the next section, we will see that such two-phonon transitions can be dominated the single-phonon 
transitions in the system, independent of the Lamb-Dickc parameter rj. 



III. MOTIONAL ROTATING WAVE APPROXIMATION 



The simplified model Hamiltonian Eq. ( |11| ) is of the form of a generalized JCM involving quadratic two-phonon 
transitions. It can be analytically solved for small squeezing parameters without external drive (fi = 0) [Q. Alter- 
natively when \v — fi| <C £l or \2v — 0| <^ fi, and g and ( C n, an approximate analytic solution can be obtained 
by an explicit diagonalization upon the elimination of the rapidly oscillating terms, i.e. the application of MRWA. In 
the interaction picture, Eq. ( pi] ) becomes 
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There exist two resonance conditions, one at Q — v where the validity condition for MRWA becomes (ryi' s /2)e 1 / 4 <C 
Vgf 1 ! 2 H], i.e. r\ <C 2e x / 4 . When the difference in trap frequencies is small we get rj <C 7/4 + (y e j2vg) 2 . We see the 
validity regime of MRWA depends on r\ for this case of single phonon JCM. Physically the trap frequency difference 



can be controlled experimentally Jl5[. When 



as for an ion, we recover the condition rj <SC 2 as discussed by 



Moya-Cessa et al. H. Under this condition, we arrive at the single phonon model 
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This Hamiltonian is the same with the Jaynes-Cummings model i/jcM 
rotating terms a<<r+,cur- are dropped under RWA. 
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Alternatively, there is a second resonance at fl — 2v. The condition for MRWA in this case is C/2 < n @, which 
becomes lv 2 + 9fg 0. Hence, the validity regime of MRWA becomes independent of ?/ in this case and can be 
realized simply by choosing larger trap frequencies. In contrast to the single phonon JCM in Eq. (|l3|), the two phonon 
squeezing model can be achieved beyond LDL and is described by 

£(2p) _ y (2p) F ( 2p ) 
n MRWA — rL hcc ' v > 

<lv( 2 P) t 

™frcc = VCL a + W «i 

y(2p) =i( t 2 o-_ + o 2 o- + ). (14) 

Equations Eq. ( |l3| ) and Eq. ( |l4| ) consist the major results of the general unitary transform introduced earlier. Their 
complete dynamics, however, can be complicated as transformation back into the original Schrodinger picture induces 
coherent mixing of different states. In particular, if the system is initially prepared in a state ip(0), the linearized 
model Hamiltonian Eq. ( [i"3| ) has an initial state ^(0) = Tfy(0), which is always a superposition of both internal states. 
In other words, even for systems initially in only one internal state such that it has a vanishing dipole moment, the 
transformed system will always have a non- vanishing dipole- moment. It is known in this case that counter rotating 
terms (CRT) neglected in making MRWA become more significant and validity regimes for MRWA more restricted 

We are now in a position to discuss validity regimes of the commonly used MRWA for the simplified model Hamil- 
tonians Eqs. (|l^) and (|l4|). We will focus on the single-phonon model in Eq. ( |l3| ) as it can be compared directly 
to existing results (under MRWA) in literature ||. Furthermore, while the single phonon JCM can be obtained for 
small restricted values of the Lamb-Dicke parameter, two-phonon model is obtained without any restriction on the 
Lamb-Dicke parameter. Thus, determining the validity regime of MRWA is more essential for the single-phonon JCM. 



IV. NUMERICAL METHOD 



In the following discussion we focus on the single-phonon model and thus ignore the unnecessary superscript (lp) 
for notational simplicity. In the numerical studies, for the approximate JCM Eq. ( |l3| ) we propagate any given initial 
states according to -ip(t) = e~ l7ii " !ot Ui(t)'tp(0) with the known propagator Uj(t) = er %w in the interaction picture [[L7| , 
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Transformed states are denoted by ip — Tip. In the Schrodinger picture, where observables are computed, the wave 
function becomes ip(t) — 'e _tWtrcc *[/ j(t)Tip(0) . For example, the mean phonon number (a^a) is found to be 

(n) =^(0) t f t L r )e mf '-*fa t af t e- mf - t C/ 7 (i)rV'(0), (16) 
which can be further simplified with the use of displacement operator property 

2 

fa)af ] = (a f a+ —)I- i|(a f - a)a x . (17) 

We note that the last term in the above Eq. involves momentum operator. Therefore, it does not commute with Hq 
and introduces fast oscillations due to the e ~ l2ut factor. We finally get, 
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We immediately see that CRT contribute to the dynamics of (n) because of multiple transitions in the original 
Hamiltonian. In addition to (n), CRT also affect other dynamical variables of interest, e.g. the Mandel-Q factor 
(((a^a) 2 ) — (a^a) 2 )/2, which characterizes phonon statistics and the impurity factor X = 1 — tr(/?p) [Q. p p is the 
reduced density operator in the phonon subsystem. 

When all CRT are included, i.e. no MRWA is made, it is no longer possible to solve dynamic propagation analytically. 
However, the transformed Hamiltonian can be diagonalized numerically in a truncated phonon Fock state basis. It is 
easy to check the accuracy of such a truncation by testing for convergence with successively larger basis states. For an 
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initial Poisson distribution of vibrational phonons, the dimension of the truncated space becomes small enough that 
efficient algorithms are readily available. To illustrate this, let us consider an initial state where the particle is in the 
internal state |e) with a motional coherent state \a), i.e. tp(0) = \e)\a). The transformed initial condition becomes 

#)) = e-*)M_^| a (19) 
In a truncated Fock space of dimension N (3> \a — ir]/2\ 2 ), the state vector is expanded as 

N 

i,{t)= X n+nA (t)\An), (20) 

n— 1, A— e : g 

with riA = 0, ./V and A = e, g respectively. The initial condition (|l9|) then becomes 

X n+nA (0) = s A ^=e-^ R ^F n ^(a - (21) 

with sa = +1, —1. The coherent state probability amplitudes in Fock space are F n (a) — exp (— \a\ 2 /2)a n / *Jn\. After 
back transformation, the original state vector evolution is then determined by 

V>(t)= A n +n K {t)\kn), (22) 

n — ljA= e,g 



with coefficients Aj for (j = 1, 2, • • • , 2N) given by 

1 N 

A n +n A (t) — —j= Yj Dn-l,m-l(sA~)[X m (t) + S A X m+N (t)}. (23) 
* m— 1 

The displacement operator matrix elements in Fock basis are known analytically in terms of Laguerre polynomials 
p9| . Therefore, the complete time evolution is obtained provided the transformed Hamiltonian is diagonalized. This 
is achieved through the standard expression 

2N 

*i=i,...w(*) = E V l3 V^e- m ^X k {Q), (24) 

where Wj are the eigenvalues of the transformed Hamiltonian and V is the diagonalizing matrix whose columns are the 
corresponding eigenvectors. In the following section we will compare MRWA results with the numerical diagonalization 
method (NDM). 



V. RESULTS AND DISCUSSIONS 



We focus in this section on the comparison of exact numerical results with results from MRWA for an initial 
motional coherent state |a) with a = \a\ exp (i/3). Moya-Cessa et al. |J have discussed a restricted form of the unitary 
transformation that led to a similar linearized model Hamiltonian earlier. Although formally independent of the 
Lamb Dicke parameter, their results were obtained under MRWA. In this first example, we demonstrate that their 
illustrative figure as presented in Ref. j| is in fact invalid at the presumed marginal LDL 77 = 0.5. In Fig. [I], we see 
that MRWA result predicts regular behaviors for several dynamical variables, while the exact result from the NDM 
at 77 = 0.5 and a — (0.5, 5) display significant differences. In general we find that with NDM, the initial state losses 
its purity faster and ends up with a larger value for the time averaged mean phonon number. On the other hand, 
MRWA results typically give larger widths of temporal fluctuations for the Mandel-Q factor, especially at earlier 
times. Perhaps most importantly, MRWA results predict super-revivals || in Q and (n) that were never observed 
with NDM. 
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Time 

FIG. 1. Comparison of the mean phonon number (n), the Mandel-Q factor, and the impurity parameter T for r\ = 0.5 and 
a = (0.5, 5). Always, the upper sub-figures are from the exact numerical diagonalization, while the lower ones are from MRWA. 
Time axis is in dimensionless from (scaled by g = rjv/2). Note the differences as compared with Ref. [3]. 

Differences of similar orders of magnitude are also found for a = (5, 0.5) as detailed in Fig. 0. We also note the 
significantly improved quantitative agreement for (n) in this case. In fact, such improved agreement with MRWA 
results always seem to occur when /3 « 0. 
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FIG. 3. The same as in Fig. 1 but for r\ = 0.1 and a — (0.5, 5). Note the super-revivals for Q. 

Figure || compares early time dynamics for the same rj = 0.1 but with a = (5, 0.5). We note that the X — value 
predicts the existence of Schrodinger cat states from both MRWA and NDM results |L8| . Also the agreement for (n) 
and Q are better than the case of a = (0.5, 5). These results show clearly that errors in MRWA are sensitive to the 
phase of the initial coherent state amplitude. 
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FIG. 4. The same as in Fig. 1 but for r] = 0.1 and a = (5, 0.5). 



To further elucidate effects of MRWA, we have carried out additional comparisons. In Fig. ||and in Fig. [|, we notice 
that MRWA results predict larger oscillations than the actual smoother behavior of (n), Q, and /. This observation 
is in contrast to our expectations that stem from the effects of RWA in the usual JCM. Let us first recall what is 
the effect of neglecting counter rotating terms in the Hjcm- For comparation, we consider an initial condition of 
the form Tip(0). The results obtained by propagating Hjcm dynamically with and without counter-rotating terms 
are given in Fig. |[ We see Hjcm with RWA predicts smoother behavior while the actual Hjcm results carry 
small nutations due to counter-rotating effects. For the ion-trap system displayed in Fig. the results are just the 
opposite. In ion-trap Hamiltonian MRWA results carry more nutations than the actual smoother behavior. The basic 
difference between the JCM and the ion-trap system is the requirement of the additional back transformation for 
the latter via T' in determining the wave function evolution. From the physical point of view, this transformation 
brings back the effect of multiplc-phonon transitions, present in the ion-trap case, after a simplified single phonon 
transition dynamics has been determined conveniently in the interaction picture. Neglecting counter-rotating terms 
in the effective single- phonon model Eq. (|13|) , however, ignores too many multiple-phonon processes at the end of the 
transformation, thus introduce more noise around the otherwise smoother collapse-revival patterns. We note that this 
(MRWA in ion-trap models) is opposite to effects of the RWA in the JCM and has been found for an initial motional 
state with almost real coherent state amplitudes, i.e. (3 » 0. 
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FIG. 5. The mean phonon number (n), the Mandel-Q factor, and the impurity parameter T with r\ = 0.1 and a = (3,0) 
for the simple single-phonon JCM Hjcm with and without RWA. In both cases no back transformation is employed but an 
initial condition of the form fip(0) has been used to compare with the multiple-phonon transition effects in Fig. ^. Within 
each sub-figure, the upper (lower) part is for J/jcm with (without) counter rotating terms. 
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FIG. 6. The same as in Fig. 5 but now for the full ion-trap Hamiltonian which includes multiple phonon transition effects. 
Within each sub-figure, the upper (lower) curve is obtained by NDM (MRWA). 



VI. CONCLUSION 

We have investigated the regimes of validity for MRWA of a trapped particle under its coherent interaction with a 
plane wave laser field. Our study is based on a family of general unitary transformations that incorporate several earlier 
transforms as special cases |j|,f^] . This general transformation facilitates the linearization of the system Hamiltonian 
including the particle's motional degrees of freedom. In model studies presented here, two distinct dynamical regimes 
appear, the single-phonon dynamics and the double-phonon dynamics. The single phonon model is standard in the 
deep LDL. The two-phonon model, on the other hand, does not impose any restrictions on the Lamb-Dicke parameter. 
Thus, addressing validity regimes of MRWA is only necessary for the single-phonon dynamics used extensively in the 
literature. For this aim, we have employed a numerical diagonalization procedure (without MRWA) and comparatively 
assessed the validity conditions of MRWA. We find that remarkable errors exist in regimes where MRWA was believed 
to be applicable. Furthermore, our simulations show that the accuracy and detailed quantitative effects of MRWA 
depends on the phase of the initial motional coherent state. When this phase is close to zero, time averaged mean 
phonon number (n) displays improved agreement with MRWA results. The agreement however is far from perfect 
even in the deep LDL. Qualitatively, for almost real initial coherent state amplitude, the counter rotating terms in 
ion-trap model lead to smoother dynamical behaviors than the results obtained under MRWA, which introduces more 
larger nutations around the actual more smooth dynamical patterns. This is in stark contrast to CRT effects in usual 
JCM problems. Physically, this effect arises because of the existence of multiple transitions in motional states. Our 
study indicates that the connection between cavity QED and ion-trap systems in JCM-type formalisms are possible 
as long as ranges of the Lamb-Dicke parameter are carefully analyzed to enforce the validity of MRWA. 
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